PASJ: Publ. Astron. Soc. Japan , 1-??, 
© 2008. Astronomical Society of Japan. 



Astrometry of Galactic Star Forming Region Sharpless 269 with VERA 
: Parallax Measurements and Constraint on Outer Rotation Curve 

Mareki Honma, 1,2 Takeshi Bushimata, 1 ' 3 Yoon Kyung Choi, 1,4 Tomoya Hirota, 1 Hiroshi Imai, 5 
Kenzaburo Iwadate, 6 Takaaki Jike, 6 Osamu Kameya, 2,6 Ryuichi Kamohara, 1 Yukitoshi Kan- ya, 1,7 
Noriyuki KAWAGUCHI, 1 ' 2 Masachika KlJlMA, 1 ' 2 Hideyuki KOBAYASHI, 1 ' 3 ' 4 ' 6 Seisuke Kuji, 6 
Tomoharu Kurayama, 1 Seiji Manabe, 2,6 Takeshi Miyaji, 1 ' 3 Takumi Nagayama, 5 
Akiharu Nakagawa, 5 Chung Sik Oh, 1 ' 4 Toshihiro Omodaka, 5 Tomoaki Oyama, 1 Satoshi Sakai, 6 
Katsuhisa Sato, 6 Tetsuo Sasao, 8,9 Katsunori M. Shibata, 1 ' 3 Motonobu Shintani, 5 
Hiroshi Suda, 4 ' 6 Yoshiaki Tamura, 2,6 Miyuki Tsushima, 5 and Kazuyoshi Yamashita 1,2 
q ' 1 Mizusawa VERA Observatory, NAOJ, Mitaka, Tokyo 181-8588 

, 2 Graduate University for Advanced Studies, Mitaka, Tokyo 181-8588 

• 3 Space V TBI Project, NAOJ, Mitaka, Tokyo 181-8588 

C*) 4 Department of Astronomy, University of Tokyo, Bunkyo, Tokyo 113-8654 

■ ''Faculty of Science, Kagoshima University, Korimoto, Kagoshima, Kagoshima 890-0065 

C/} \ 6 Mizusawa VERA observatory, NAOJ, Mizusawa, 023-0861, JAPAN 

7 Department of Astronomy, Yonsei University, Seoul, 120-749, Republic of Korea 
s Ajou University, Suwon 44^-^49, Republic of Korea 
9 Korean VLBI Network, KASI, Seoul, 120-749, Republic of Korea 
honmamr@cc.nao.ac.jp 



oo 
o 



- I - < 

>< 



Q \ (Received 2007 March 9; accepted 2007 May 20) 

5— I ' 

■ Abstract 

CZ2 

i We have performed high-precision astrometry of H2O maser sources in Galactic star forming region 

Sharpless 269 (S269) with VERA. We have successfully detected a trigonometric parallax of 189 ± 8 /Ltas, 
corresponding to the source distance of 5.28±g;^ kpc. This is the smallest parallax ever measured, and 
the first one detected beyond 5 kpc. The source distance as well as proper motions are used to constrain 
the outer rotation curve of the Galaxy, demonstrating that the difference of rotation velocities at the Sun 
and at S269 (which is 13.1 kpc away from the Galaxy's center) is less than 3 %. This gives the strongest 
constraint on the flatness of the outer rotation curve and provides a direct confirmation on the existence 
of large amount of dark matter in the Galaxy's outer disk. 

Key words: ISM:star forming regions — ISM:individual(Sharplcss 269) — mascrs (H2O) — VERA 

1. Introduction molecular clouds, and 2) in most cases it has been diffi- 

cult to measure proper motions of Galactic objects, and 
Rotation curves, which are the plots of rotation veloci- thus only radial velocity (out of the three components of 



ties as functions of distance from galaxy centers, arc im- spatial velocity) could be used to constrain the Galactic 

■ portant tools to study the mass distributions in disk galax- rotation. Therefore, outer rotation curve and hence dark 

ies. Assuming the centrifugal force balances with gravity, matter distribution in the Galaxy's disk is still highly un- 

one can determine the mass distribution in a galaxy us- certain, although it is widely believed that the Galaxy's 

ing the rotation curve. Observations of rotation curves outer rotation curve is flat just like those of extra-galaxies, 

in external galaxies revealed that they are basically flat In order to determine a precise rotation curve of the 

within (and often beyond) optical disks of spiral galax- Galaxy, high-precision astrometry is essential. By mea- 

ies (e.g., Rubin ct al. 1980, Rubin 1983, Sofue & Rubin suring the accurate position of a star, one can deter- 

2001). Flat rotation curves observed in external galax- mine a trigonometric parallax tt and thus source dis- 

ies provided strong evidences for the dark matter in the tance as D = I/71". In addition, high-precision astrom- 

galaxies' outer regions (e.g., Kent 1986, 1987). etry also allows us to determine source proper motions 

In contrast, the rotation curve of the Milky Way Galaxy (source motions projected onto the sky plane), and thus 

remains highly uncertain particularly in the outer region, the 3-dimensional space velocity of the source. Yet accu- 

although a large number of efforts have been made to de- rate astrometric measurements have been done only for 

termine it (e.g, Clemens 1985, Mcrrificld 1992, Brand & sources that are located fairly close to the Sun compared 

Blitz 1993, Honma & Sofue 1997, and references therein), to the Galaxy's size. For instance, the HIPPARCOS mis- 

The reason for this is two fold: 1) it has been difficult to sion (Perryman et al. 1997), the most modern satellite 

measure accurate distances of the Galactic objects which dedicated to optical astrometry, has reached distances of 

are used to trace the rotation curve, such as OB stars and ~300 pc by parallax measurements, which is much smaller 
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than the size of the Galaxy (e.g., ~15 kpc for the radius 
of the Galaxy's disk). Hence, at present astrometry of the 
Galaxy still remains an unexplored issue. 

Over the next decades, there will be new missions 
for astrometry of the Galaxy that aim at 10 /Ltarcsec 
accuracy (e.g., SIM 1 , GALA 2 , JASMINE 3 ). These are 
satellite- type missions with optical/infra-red telescopes 
orbiting the Earth, where the observations can be made 
free from the disturbance of the atmosphere. At radio 
wavelengths, high precision astrometry have been per- 
formed with ground-based VLBI (Very Long Baseline 
Intcrfcromctry) , a radio interferometer with telescopes 
separated by ~ 1000 to ~ 10000 km. The advantage of 
VLBI is that it provides the highest angular resolution 
among the existing telescopes at any wavelength. While 
normal VLBI observations directly suffer from the fluc- 
tuation of atmosphere (mainly due to the water vapor 
content in the troposphere), a noble way, which is called 
'phase-referencing', has been developed to cancel out the 
tropospheric fluctuations. In phase-referencing observa- 
tions, a few sources (one target and one or more refer- 
ence sources) are observed at nearly same time by rapidly 
switching the telescope, and then relative positions of the 
target with respect to the reference sources can be mea- 
sured after correcting for the influence of troposphere. In 
fact, recently the distance of the Galactic star forming re- 
gion W30H was determined with VLB A (Xu ct al. 2006, 
Hachisuka et al. 2006) to solve the long-standing ambigu- 
ity of the Perseus arm distance by measuring the source 
distance at 2 kpc. It has also been shown that the maser 
emissions from late-type stars (such as Mira variables) 
can be used to perform kpc-scale astrometry with VLBI 
(Kurayama et al. 2005). These results strongly demon- 
strate the promising future of ground-based VLBI for the 
Galaxy-scale astrometry. 

VERA (VLBI Exploration of Radio Astrometry) is a 
new Japanese VLBI array dedicated to VLBI astrometry 
(Honma et al. 2000, Kobayashi et al. 2003). A unique 
feature of the VERA telescopes is the dual-beam receiv- 
ing system, which allows us to observe simultaneously a 
target and a reference source within 2.2 degree. This dual- 
beam system enables to cancel out tropospheric fluctua- 
tions much more effectively than switching observations 
with normal, single-beam telescopes (Honma et al. 2003). 
With a target astrometric accuracy of 10 ^arcsec, VERA 
aims to precisely locate hundreds of maser sources in the 
Galaxy and to explore its 3D structure and dynamics. The 
VERA array was constructed by 2002, and regular obser- 
vations have been started in the fall of 2003. Here we 
present initial results of high precision astrometry with 
VERA toward the Galactic star forming regions S269, 
which is located in the outer regions (/ = 196.45°). We 
report on the determination of the smallest parallax in 
human history. Our results provide the strong constraint 
on the flatness of rotation curve in the outer Galaxy. 

1 SIM: http://planctqucst.jpl.nasa.gov/SIM/sim_indcx.cfm 

2 GAIA: http://sci.esa.int/science-e/www/area/index.cfm? 
farcaid=26 

3 JASMINE: http://www.jasmine-galaxy.org/indcx.html 



2. Observations and Reductions 

We observed H2O masers in Galactic star forming re- 
gion Sharplcss 269 (S269) with VERA since November of 
2004 and here we present the data of 6 epochs that were 
obtained with the full 4-station array (Mizusawa, Iriki, 
Ogasawara, and Ishigaki-jima) under relatively good con- 
ditions. The epochs are day of year (DOY) 323 in 2004, 
DOY 026, 073, 134, 266 and 326 in 2005 (Nov. 18 in 
2004, Jan. 26, Mar. 14, May. 14, Sep. 23 and Nov. 21 
in 2005), spanning ~1 year. At each epoch, H2O 6i6-523 
maser line at a rest frequency of 22.235080 GHz in S269 
and a position reference source J0613+1306 were simulta- 
neously observed in dual-beam mode for nearly 9 hours. 
Typical on-source integration time was 5 hours for both 
the target maser and the reference. The reference source, 
J0613+1306, is one of the ICRF sources (Ma ct al. 1998) 
with a correlated flux density of ^300 mJy The sepa- 
ration angle between the maser and reference sources is 
0.73 degree. Left-hand circular polarization signals were 
received for both S269 and J0613+1306, and digitally 
recorded onto magnetic tapes with the VERA-terminal 
system at the total data rate of 1024 Mbps. With 2-bit 
quantization, this data rate provides a total bandwidth of 
256 MHz. The signals from the two sources are filtered 
with VERA Digital Filter (Iguchi ct al. 2005) to obtain 
1 IF (Intermediate Frequency) channel of 16 MHz for the 
S269 maser line and 15 IF channels of 16 MHz (240 MHz 
in total) for J0613+1306. Correlation processings were 
made with the Mitaka FX correlator. For the reference, 
which is a continuum source, the spectral resolution was 
64 points per each 16 MHz channel, which corresponds to 
a velocity resolution of 3.4 km s . For the maser source, 
frequency and velocity resolutions were 15.625 kHz and 
0.21 km s , respectively. 

Since the correlator's a priori delay model is not ac- 
curate enough for high precision astrometry, recalcula- 
tions of precise delay were made after the correlation, 
and correlated visibilities were corrected for the differ- 
ence between the first (rather crude) a priori model and 
second (more accurate) delay model. The delay recalcula- 
tion code is based on the geodynamics models described in 
IERS convention 1996 (McCarthy 1996), and earth orien- 
tation parameters (EOP) were taken from IERS bulletin 
B final values 4 , which currently provides the best esti- 
mates. Also, ionospheric delays were corrected based on 
the global ionosphere map (GIM), which arc produced by 
the University of Bern every day 5 . 

In the data analysis of each epoch, at first fringes were 
searched for the reference source J0613+130. With 240 
MHz bandwidth and typical system noise temperature of 
200 K, the reference source, having a flux of ~300 mJy, 
was easily detected within 1 minute integration, which is 
much shorter than typical coherence time at 22 GHz (2 
to 3 minutes). Since its position is accurately known, the 
fringe parameters of J0613+1306 are also used to calibrate 

4 http://hpicrs.obspm.fr/eop-pc/ 

5 http:/ /www. aiub.unibe.ch/ionospherc. html 
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Fig. 1. (a, top) Auto-correlation spectrum of S 269 H2O 
maser emissions taken on DOY 073 in 2005. (b, bottom) 
maser spot distribution obtained with VERA on DOY 073 in 
2005. The coordinates are with respect to the tracking center 
position, which is (06hl4m37.08s, +13d49'36.7") in J2000. 

clock offset parameters (such as delay and delay-rate off- 
set). The phase solutions for J0613+1306 were converted 
into the phase at observed maser frequency, and applied to 
the visibilities of S269 together with the dual-beam phase 
calibration data. This dual-beam phase calibration data 
was taken real-time during the observations and are based 
on the correlation of artificial noise sources injected into 
two beams at each station (Kawaguchi ct al. 2000). After 
those calibrations, the visibilities of the S269 masers were 
Fourier transformed to synthesize images, and the posi- 
tions of the brightness peaks were determined with respect 
to the reference spot. In some epochs (especially in sum- 
mer) , the qualities of phase- referenced maps were not high 
due to residuals in troposphcric delay. To calibrate them, 
residual zenith delays were estimated as a constant off- 
set that maximizes the coherence of the phase-referenced 
map. Typical residuals of zenith delay are 1 to 5 cm, 
but in the worst case (during the summer at Ishigaki-jima 
station) it was as large as 20 cm. 

3. Results 

The total-power spectrum of S269 taken on DOY 073 
in 2005 is shown in figure 1(a). Basically it consists of a 



single feature at Vlsr 6 of ~ 19.6 km s _1 with a peak inten- 
sity of 480 Jy, being consistent with previous single-dish 
monitoring study (Lekht 2000). This main feature was 
always bright and observable for all the epochs presented 
here. Figure 1(b) is the maser spot map of the main fea- 
ture around Vlsr of ~ 19.6 km s" 1 (for DOY 2006/073). 
Six maser spots were detected in the velocity range from 
19.0 to 20.1 km s _1 , and these maser spots are aligned in 
the east-west direction on a scale of 0.4 mas. It is remark- 
able that the thickness of the feature (spots distribution 
in the north-south direction) is ~ 50 fias, 10 times smaller 
than the width in the east-west direction. The maser dis- 
tribution also shows a velocity gradient from the east to 
the west. This kind of structure is rather unusual for 
H2O masers in star forming regions (which mostly show 
bipolar structure with discrete blue/red-shifted clusters 
of maser spots), but instead, similar to those of CH3OH 
(methanol) masers at 6.7 GHz and 12.2 GHz in terms 
of the spot distribution as well as the velocity width 
(e.g., Minner et al. 2000). Note that positions in fig- 
ure 1(b) are the residuals to the tracking center positions 
of the maser and reference sources, which are taken to 
be (a, <5) = (06hl4m37.08s, +13d49'36.7" ) for S269, and 
(06hl3m57.692764s, + 13d06'45.40116") for J0613+1306, 
both in the J2000 coordinates. Thus, the absolute po- 
sition of the brightest 19.6 km s" 1 spot at DOY 073 of 
2005 is obtained as (06hl4m37.07933s, +13d49'36.6945j 
with an uncertainty of 1 mas, which mainly comes from 
the uncertainty of absolute position of the reference source 
J0613+1306. The absolute position of the maser feature 
shown in figure 1(b) agrees well with the position of S269 
IRS2w, which is the most luminous infrared source in the 
S269 regions (Jiang et al. 2003). 

Figures 2 show the positional variations of the bright- 
est maser spot in X (= cos<5Aa, east- west offset) and Y 
(= A<5, north-south offset) directions for a monitoring span 
of 1 yr. As clearly seen from the plot (especially in X 
direction), the positions show systematic sinusoidal mod- 
ulation with a period of 1 yr. The phase of observed sinu- 
soidal curve (i.e., the peak date) perfectly matches with 
that of the expected parallax curve for S269, ensuring that 
the modulation certainly originates from the parallax of 
S269. In figure 2, error bars were estimated as the stan- 
dard deviation from the best fit with parallax plus linear 
proper motions. This error estimate was made because 
it is difficult to predict the observational error in VLBI 
astrometry: the error depends on many factors such as 
residual phase in phase-referencing, error in zenith delay 
of troposphere and ionosphere, and error in calibration of 
instrumental offset and so on, and is hardly predictable. 
The estimated error bars are 25 /xas for X and 75 /xas for 
Y. We note that error in Y is three times larger than 
that in X. This can be explained if one assume that ma- 
jority of the error comes from the uncertainty in the tro- 
posphcric zenith delay, as the troposphcric delay changes 



the LSR velocity described here is in traditional definition using 
the Solar motion of (U, V, W)=(10, 15.4, 7.8) km s -1 (Kerr & 
Lyndcn-Bcll 1986) for comparisons with previous observations. 
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Table 1. The best fit values of parallax n and proper motions [ix and fiy for the three brightest spots. The last row 
gives the values obtained by weighted mean. Note that positions X and Y are those at the first epoch (DOY 323 in 
2004) with respect to the tracking center position of S269, which was taken to be (06hl4m37.08s, +13d49'36.7" ) in J2000. 

X (mas) Y (mas) Vlsr t (mas) fix (mas/yr) /iy (mas/yr) 

^8l83 1^571 1^8 0.208±0.030 -0.425±0.032 -0.123±0.076 

-8.365 -5.591 19.6 0.199±0.012 -0.388±0.014 -0.118±0.071 

-8.427 -5.602 19.3 0.176±0.012 -0.457±0.015 -0.123±0.072 

mcaH 0.189±0.008 -0.422±0.010 -0.121±0.042 
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Fig. 2. Positional variations of the brightest maser spot at 
Vlsr = 19-6 km — 1 . Top is for X (east-west) direction and 
bottom is for Y (north-south) direction. The time origin is 
the first epoch of our observations, which is DOY 323 in 2004. 
Solid curves are best-fit results with parallax and proper mo- 
tions, and dashed lines show the best-fit proper motions. 

apparent elevations of the source. Detailed consideration 
on the error is discussed in the next section. For our par- 
allax measurements, in this paper we consider only the 
X component because the V-dircction error is large and 
also because S269 is near the ecliptic and the parallax el- 
lipse is highly elongated in the X direction, making the 
contribution of Y to the parallax determination smaller. 
Here we use the brightest three spots at the radial veloc- 



ity from 19.4 to 19.8 km s" 1 including the brightest one 
at 19.6 km s _1 to ensure the high signal-to-noise ratio: 
as clearly seen from figure 1, the H2O maser spectrum is 
sharply peaked at 19.6 km s _1 , and maser intensity be- 
comes weak at off-peak radial velocities. Least-squares fits 
were made to positions of the three maser spot, with par- 
allax 7r as well as proper motions fix (east-west direction) 
and [iy (north-south direction). Table 1 summarizes the 
best fit results. As seen in table 1, independent analyses 
for the three spots give results consistent with each other, 
from 176 fias to 208 /ias. To obtain the best estimate of 
the parallax w, we took a weighted mean of the parallaxes, 
yielding the parallax of S 269 as 189 ± 8 /ias, where the er- 
ror bar is estimated as tr 2 = l/'C(l/ f )• This corresponds 
to a source distance of 5. 28^22 kpc. This is the smallest 
parallax ever measured to date, demonstrating the high 
capability of VERA to perform Galactic-scale astrometry. 
The distance to S269 is found to be slightly larger than 
previous estimates, which claimed a distance of ~4 kpc 
(Moffat et al. 1979). 

From the fitting results in table 1, the weighted means of 
the proper motions are obtained as (fix = —0.422 ± 0.010, 
fiy = —0.121 ± 0.042) mas yr _1 , respectively. To convert 
these observed (heliocentric) proper motions to the ones 
with respect to LSR, we use the solar motion based on 
the HIPPARCOS satellite data (Dehnen & Binney 1998), 
which is (U, V, W)=(10.0, 5.25, 7.17) km s" 1 . Using the 
galactic coordinates of S269 (I, 6)=(196°.45, -1°.96) and 
the Galactic plane's position angle of 151.51° there, one 
can calculate the proper motion projected to the direction 
of I and b as (fi u fi b )={-0. 184 ± 0.032, -0.149 ± 0.029) 
mas yr _1 . Given the source distance of 5.28 kpc, these 
proper motions correspond to the velocity vector of (vi, 
u 6 )=(-4.60 ± 0.81, -3.72 ±0.72) km s" 1 , respectively. 
These velocity components are remarkably small com- 
pared to the rotation speed of the Galaxy, which is an 
order of ~ 200 km s _1 . Given that S269 is located in 
the anti-center region, the small value of vi indicates that 
the Galactic rotation velocities at the Sun and at S269 
are close to each other and proper motions were cancelled 
out in our relative proper motion measurements. Detailed 
discussion on the Galactic rotation velocity will be done 
in next section. 

4. Discussion 

4--1- Sources of Astrometric Error 

As described in the previous section, the astrometric 
errors estimated from the fitting deviations are 25 /ias for 
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Fig. 3. Plot of radial velocity v r and tangential velocity iij with respect to LSR as a function of the rotation velocity ratio 8 = 0/0o. 
Observed value of vi(=— 5 ± 5) is consistent with the rotation velocity ratio of 8 = 1.00 ±0.03, showing the rotation velocity is the 
same at the position of S269 and at the Sun. Note that v r (=20 ± 5) can be explained consistently with 8 = 1.00 ± 0.03. 



X and 75 [i&s for Y. This can be explained if the dom- 
inant error source is the uncertainty in the tropospheric 
zenith delay. For instance, if we take a typical uncer- 
tainty of 3 cm for the the tropospheric zenith delay, then 
it causes an path length difference of 0.4 mm (=30 mm 
x 0.7 deg / 57.3 deg/rad, where 0.7 deg is the separation 
angle between S269 and the reference source) between the 
two sources. This uncertainty in the path length differ- 
ence roughly corresponds to 40 [i&s (=0.4 mm / 2.3xl0 9 
mm, where 2.3 xlO 9 mm is the maximum baseline length 
of the VERA array). In practice, observations are not 
done toward the zenith, but sources are at lower eleva- 
tion angle (EL). As our observations were usually done 
for elevation angle larger than 20°, the effect of zenith 
delay error is multiplied by a factor of 1 (corresponding 
to EL = 90°) to -3 (EL = 20°) depending on the source 
elevation. If this factor is taken into account, one can ex- 
pect an astrometric error of 40 to 120 /ias depending the 
EL distribution, and the astrometric error in Y obtained 
above (75 /zas) is certainly in this possible range. On the 
other hand, the astrometric error in the X direction can 
be suppressed for two reasons: first, the source pair con- 
sidered here has a smaller separation in X direction than 
in Y direction, and second, the observational track in each 
epoch is roughly symmetric with respect to the meridian 
transit (i.e., each epoch has nearly same track before and 
after the transit), and this symmetry can help to reduce 
the astrometric error in the X direction caused by the 



tropospheric zenith delay offset. 

Other possible sources of astrometric error are those in 
station positions, delay model, and ionosphere. Currently, 
VERA station positions arc determined with an accuracy 
of ~ 3 mm based on geodetic observations at S/X (2/8 
GHz) bands carried out every two weeks. Also, our delay 
calculation code is compared with CALC 7 developed by 
the NASA/GSFC VLBI group, which is the international 
standard of delay model, and it turned out that the differ- 
ence between the two codes is less than 2 mm. Therefore, 
these errors are smaller than that of zenith delay by an or- 
der of magnitude. Regarding the ionosphere, its contribu- 
tion is small at 22 GHz, and corrections with GPS-based 
Global Ionosphere Map (GIM), provided by University 
of Bern, are precise enough for 10 fias astrometry. We 
note that the trend of the larger error in Y direction was 
also found in other observations of VERA. Therefore, the 
larger error in Y is not a special phenomenon of only the 
S269 observations but rather common to all of VERA ob- 
servations, and is most likely to originate from the tropo- 
spheric zenith delay error. 

4-2. Constraint on Galactic Rotation 

Here we use the proper motions and parallax obtained 
in the previous section to constrain the Galactic rota- 
tion velocity at the position of S269. First, the small 
proper motion perpendicular to the Galactic plane (vb = 

7 http:/ /gemini. gsfc.nasa.gov/solve/ 
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Fig. 4. Rotation curve of the Milky Way Galaxy obtained in previous studies together with our result on S269 using Galactic 
constants of Ro=8 kpc and So =200 km . Dashed line is the flat rotation curve with 0=200 km , shadowed area shows 
possible range of outer rotation curves considered in previous studies (Honma & Sofue 1997). Points at R< 8 kpc are inner rotation 
curves determined from the tangential velocities of Galactic HI gases (Honma & Sofue 1997), with a smoothed fit (thick curve). 
Dotted curve is the rotation curve for an exponential disk, corresponding to a constant mass-to-light ratio disk without dark matter. 
Discrepancy between the observed point for S269 and the exponential disk is evident, demonstrating the existence of large amount 
of dark matter in the outer region of the Galaxy. 



—3.72 ±0.72 km s _1 ) indicates that S269 basically par- 
takes in galactic rotation and also that the H2O maser 
proper motions truly reflect the systemic motion of S269. 
Radial velocity measurements supports the latter idea: 
while the maser emission is peaked at f9.6 km , the 
peak velocity of HII region observed with SII lines is 16.5 
km s , and the systemic velocity of the associated molec- 
ular cloud observed in CO is 17.7 km s _1 (Godbout et 
al. 1997), which agrees well with the maser radial veloc- 
ity within 3 km s . Therefore, from the radial velocity 
as well as proper motions perpendicular to the Galactic 
plane, one can expect that the peculiar velocity of the 
maser source with respect to pure Galactic rotation is as 
small as ^5 km s _1 . 

From this fact one can expect that the proper motion 
in I direction (vi) basically reflects the difference of the 
galactic rotation velocities at the position of S269 and at 
the Sun, and also that one can constrain the outer rotation 
velocity at the position of S269. In fact, S269 is located 
near the galactic anti-center region (I = 196°. 45), and thus 
the lack of a large proper motion in I direction indicates 
that the galactic rotation speed there is close to that at 
the Sun. If a source is on perfect Galactic rotation, radial 
and tangential velocities with respect to LSR observers 
can be written as, 



«*=(§- 1*) ifacoBl-llJ, (2) 

For the component perpendicular to the Galaxy plane, 
obviously Vf, = 0. These equations relate the observed ve- 
locities to the rotation velocity at the source (9) through 
the Galactic constants i?o and Qq. Figure 3 shows the plot 
of expected v r and vi for S269 as a function of 8 = 0/©o 
(rather than for seeing the difference of the rotation 
velocities at S269 and the Sun). Here the Galactic con- 
stant i?o is assumed to be 8.0(±0.5) kpc (Reid 1993), and 
three cases of ©o, 180, 200, and 220 km s" 1 arc consid- 
ered as recent determinations of Oo still vary substantially 
(Oiling & Merrifield 1998, Miyamoto & Zhu 1998, Reid 
& Brunthalcr 2004). As seen in figure 3, the tangential 
velocity vi has a strong dependence on 9, and thus ob- 
served vi can be used to constrain the rotation velocity. 
Considering that the velocity component perpendicular to 
the Galactic plane is 3.7 km s _1 and that the radial veloc- 
ity differences between maser and other lines (such as CO 
and SII) are also ^3 km s _1 , here wc can safely assume 
the possible range of the true tangential velocity as u;=-5 
± 5 km s _1 and that of the true radial velocity as v r =20 
± 5 km s -1 . 8 From figure 3, this tangential velocity range 



Note that the radial velocity considered here is slightly differ- 
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is obtained if 9 lies between 0.97 and 1.03. Therefore, at 
the position of S269, 9 should be 1.00 ±0.03, i.e., the ro- 
tation velocity at the S269 (O) must be the same with Oo 
within 3% level. This is the strongest constraint of the 
rotation velocity in the outer galaxy ever obtained. Note 
that the same argument is possible based on the radial 
velocity v r , but the constraint is not as strong as that 
from the tangential velocity, because the gradient of v r (9) 
plot is much shallower. However, we note that the radial 
velocity v r — 20 ±5 km s _1 can be consistently explained 
by the rotation velocity ratio of 9 ~ 1. 

In previous works, the Galactic rotation curve has an 
uncertainty up to 100 km s" 1 in the outer region if one 
includes the strong dependence on Galactic constants &q 
(Honma & Sofue 1997). This situation is summarized in 
figure 4, showing the area of uncertainty in previous stud- 
ies. The point for S269 determined in this study is also 
shown in figure 4. The coincidence of rotation velocities at 
the Sun and S269 simply indicates that the rotation curve 
there is basically flat, as was known for rotation curves 
of other spiral galaxies (Rubin ct al. 1980, Rubin 1983, 
Sofue & Rubin 2001). In disk galaxies like the Galaxy, 
optical surface brightness obeys an exponential law, i.e., 
I(r) = ioexp(— r/h), where r is the radius and h is the disk 
scale length. Assuming that surface brightness traces the 
mass density (i.e., constant Mass-to-Light ratio), one can 
calculate the rotation curve of optical disk assuming no 
dark matter (Freeman 1970). In figure 4, we also showed 
such a rotation curve for the Galaxy's disk, assuming a 
disk scale length of h = 3 kpc and a maximum rotation 
velocity of 200 km _1 . As seen from figure 4, such a rota- 
tion curve without dark matter was not completely ruled 
out previously. However, the rotation velocity measure- 
ment for S269 evidently shows the discrepancy between 
the observed rotation curve and that expected from the 
optical disk without dark matter, providing a strong con- 
firmation on the existence of dark matter in the Galaxy's 
outer region. At the position of S269, the rotation curve 
of the exponential disk gives V cxp — 168 km s^ 1 while the 
astromctric measurement of S269 provides V h s = 200 km 
s . In case of a spherical mass distribution, the enclosed 
mass M r within a radius r is proportional to V 2 , as M r can 
be written as Mr k rV 2 /G. The values of V exp and V b s 
obtained above give (V eX p/V h s ) 2 = 0.70, and thus within 
the position of S269 at least ~ 30 % of the enclosed mass 
must be composed of dark matter. 

Although this study presents astromctric measurements 
for only one source, the present study demonstrates the 
high capability of VERA in studying the Galactic rotation 
based on Galaxy-scale astrometry. During next decade, 
VERA continues astrometric observations of nearly one 
thousands Galactic maser sources and will provide an ac- 
curate rotation curve over the whole Galaxy's disk as well 
as an accurate description of the dark matter distribution 
of the Galaxy. 

ent from the traditionally-defined Vlsr. since the Solar motion 
considered here is the one recently obtained from HIPPARCOS 
data (Dehnen & Binney 1998). However, the difference in Vlsr 
and v r is not significant, being 3.3 km _1 for S269. 
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